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ABSTRACT 

The emission mechanism and the origin and structure of magnetic fields in gamma-ray burst (GRB) jets 
are among the most important open questions concerning the nature of the central engine of GRBs. In spite 
of extensive observational efforts, these questions remain to be answered and are difficult or even impossible 
to infer with the spectral and lightcurve information currently collected. Polarization measurements will lead 
to unambiguous answers to several of these questions. Recent developments in X-ray and 7-ray polarime- 
try techniques have demonstrated a significant increase in sensitivity enabling several new mission concepts, 
e.g. POET (Polarimeters for Energetic Transients), providing wide field of view and broadband polarimetry 
measurements. If launched, missions of this kind would finally provide definitive measurements of GRB polar- 
izations. We perform Monte Carlo simulations to derive the distribution of GRB polarizations in three emission 
models; the synchrotron model with a globally ordered magnetic field (SO model), the synchrotron model with 
a small-scale random magnetic field (SR model), and the Compton drag model (CD model). The results show 
that POET, or other polarimeters with similar capabilities, can constrain the GRB emission models by using 
the statistical properties of GRB polarizations. In particular, the ratio of the number of GRBs for which the 
polarization degrees can be measured to the number of GRBs that are detected (Nm/Nd) and the distributions 
of the polarization degrees (11) can be used as the criteria. If Nm/Nd > 30% and 11 is clustered between 0.2 
and 0.7, the SO model will be favored. If instead Nm/Nd < 15%, then the SR or CD model will be favored. If 
several events with 11 > 0.8 are observed, then the CD model will be favored. 

Subject headings: gamma rays: bursts — magnetic fields — polarization — radiation mechanisms: non-thermal 



1. EVITRODUCTION 

Gamma-ray bursts (GRBs) are brief, intense flashes of 7- 
rays originating at cosmological distances, and they are the 
most luminous objects in the universe. They also have broad- 
band afterglows long-lasting after the 7-ray radiation has 
ceased. It has been established that the bursts and after- 
glows are emitted fr om outflows movin g towards us at highly 
relativistic speeds dTavlor et alj |2004|) . and at least some 
GRBs are associated with the collaps e of massive stars (e.g., 
iHiorfli et al.ll2003l; IStanek et al.ll2003h . Observations suggest 
that the burst is produced by internal dissipation within the 
relativistic jet that is launched from the center of the explo- 
sion, and the afterglow is the synchrotron emission of elec- 
trons accelerated in a collisionless shock driven by the inter- 
action of the jet with the surrounding medium (for recent re- 
views, Piran 2005: .Meszaros.2006: .Zhang 2007.) . 
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In spite of extensive observational and theoretical efforts, 
several key questions concerning the nature of the central en- 
gines of the relativistic jets and the jets themselves remain 
poorly understood. In fact, some of these questions are very 
difficult or even impossible to answer with the spectral and 
lightcurve information currently collected. On the other hand, 
polarization information, if retrieved, would lead to unam- 
biguous answers to these questions. In particular, polarimetric 
observations of GRBs can address the following: 

Magnetic composition of GRB jets - It is highly speculated 
that strong magnetic fields are generated at the GRB central 
engine, and may play an essential role in the launch of the 
relativistic jets. However, it is unclear whether the burst emis- 
sion region is penetrated by a globally structured, dynami- 
cally important magnetic field, and whether the b urst is due to 
shock dissipation or magnetic reconnection (e.g.,ISpruit et alj 



1200 U IZhang & Mesz arosll2002l: iLyutikov et al.y2003l)~ 

Emission mechanisms of the bursts - The leading model 
for the emission mechanism of the prompt burst emission is 
synchrotron emission from relativistic electrons in a globally 
ordered magnetic field carried from the central engine, or ran- 
dom magnetic fields generated in-situ in the shock dissipa- 
tion region (Rees & Meszaros 1994). Other sug gestions in- 
clude Compton drag of arnbient s oft photons (iShaviv & Dad 
1995t lEichler & LevinsonI 120031: iLevinson & Eichled 120041: 
La zzati et alj 20041). synchro tron self-Compton emission 
(Pana itescu & Meszarosl20 0u), and the combination of a ther- 
mal component from the ph otosphere a nd a non-thermal com- 
ponent (e.g., synchrotr on) (fRvde et al.ir2006l: [Thompson et alj 
l2007tllokaetani2007h . 

Geometric structure of GRB jets - Although it is generally 
believed that GRB outflows are collimated, the distribution of 
the jet opening angles, the observer's viewing direction, and 
whether there are small-scale structures within the global jet 
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are not well understoo d dZhang et alj l2004t lYamazaki et alj 
12001: iToma eTaPIIOOl . 

To date, robust positive detections of GRB polarization 
have been made only in the optical band in the afterglow 
phase. Varying linear polarizations have been observed in 
several optical afterglows several hours after the burst trigger, 
with a level of ~ 1-3%, which is consistent with the syn- 
chr otron emissio n mechanism of GRB afterglow (for reviews, 
see ICo vino et al . 2004; Lazzati 2006). An upper Hmit (< 8%) 
has been obtaine d for the early (t ~ 2 00 s) optical afterglow 
of GRB 060418 jMundell et al.|[2007l) . Also for radio after- 
glows, we have several upper limits for the polarization degree 
jTavlor et al.l2005l:lGranot & Tavlod2005 ') (for some implica- 
tions. see lToma etal .I2008h . As for the prompt burst emission, 
strong linear polarization of the 7-ray emission at a level of 
n = 80 ± 20% was clai med for GRB 021 206 based on an anal- 
ysis of RHESSI data jCoburn & Boggsll2003h . although this 
claim remains contr oversial because of large systema tic un- 
certainties (RutledgeA^illOOl iMggeZSiilllQQl). Sev- 
eral other reports of high levels of polarization in the prompt 
burst emissi on are a lso statistic ally inconc lusive (Wi llis et al.l 
|2005; KalemcietaljL2007; Mc Glvnn et al.ll2007l) . 

Recently, more sensitive observational techniques for X- 
ray and 7-ray polarimetry have been developed, and there 
are several polarimeter mission concepts. These include Po- 
larimeters for Ener getic Transients (POET, Hill et al. 2008; 
iBloser et al.l l2008|). P olarimeter of G amma-ray Observe r 



(PoGO. iMizuno et all [BoS), POLAR dProduit et al.l 



Advanced Compton Telescope (ACT. iBoggs et alT 



2005) 



2006), 



Gravity and Extreme Magnetism (G/sM^. Ijahoda e t al.''2007^ 
XPOL (Costa etal. 2007), Gamma-Ray Burst Inv estigation 
via P olarimetry and Spectroscopy (GRIPS, iGreiner et alj 
l2008h . and so on. 

Several of these missions, if launched, would provide 
definitive detections of the burst polarizations and enable us 
to discuss the statistical properties of the polarization degrees 
and polarization spectra. Although there are several polarime- 
try mission concepts described in the literature, POET is the 
only one to date that incorporates a broadband capability for 
measuring the prompt emission from GRBs, and for this rea- 
son it provides a good case study for our simulations. POET 
will make measurements with two different polarimeters, both 
with wide fields of view. The Gamma-Ray Polarimeter Exper- 
iment (GRAPE; 60-500 keV) and the Low Energy Polarimeter 
(LEP; 2-15 keV) provide a broad energy range for the obser- 
vations. Suborbital versions of both POET instruments are 
currently being prepared for flight within the next few years. 
GRAPE will fly on a sub-orbital balloon in 2011, and the 
Gamma-Ray Burst Polarimeter (GRBP, a smaller version of 
LEP) will fly on a sounding rocket. 

Theoretically, it has been shown that similarly high lev- 
els of linear polarization can be obtained in several GRB 
prompt emission models; the synchrotron model with a 
globally ordered magnetic field, the sync hrotron model 
with a small- scale random magnetic field (iGranotI 120031 : 
iLvutikov etai] 12003; Nakaretal. 2003) , and the Compton 
drag m odel (Lazzati et al. 2004; Eichler & Levinson 2003; 
iLevinson & Eichler 2004: Shaviv & Dar 1995). Thus the de- 
tections of GRB prompt emission polarization would support 
these three models. In this paper, we show that these mod- 
els can be distinguished by their statistical properties of ob- 
served polarizations. We performed detailed calculations of 
the distribution of polarization degrees by including realistic 
spectra of GRB prompt emission and assuming realistic dis- 



tributions of the physical parameters of GRB jets, and show 
that POET, or other polarimeters with similar capabilities, can 
constrain the GRB emission models. We use the limits of 
POET for GRB detection and polarization measurements as 
realistic and fiducial limits. This paper is organized as fol- 
lows. We first introduce the POET mission concept in §|2l In 
§ [3] we summarize the properties of the observed linear po- 
larization from uniform jets within the three emission models. 
Based on these models, we perform Monte Carlo simulations 
of observed linear polarizations and show how the statistical 
properties of observed polarization may constrain GRB emis- 
sion mechanisms in §|4] A summary and discussion are given 
in §111 

2. PROPERTIES OF POET SATELLITE 

POET (Polarimeters for Energetic Transients) is a Small 
Explorer (SMEX) mission concept, that will provide highly 
sensitive polarimetric observations of GRBs and can also 
make polarimetry measurements of solar flares, pulsars, soft 
gamma-ray repeaters, and slow transients. The payload con- 
sists of two wide field of view (FoV) instruments: a Low En- 
ergy Polarimeter (LEP) capable of polarization measurements 
in the 2-15 keV energy range and a high energy polarime- 
ter (Gamma-Ray Polarimeter Experiment; GRAPE) that will 
measure polarization in the 60-500 keV energy range. POET 
can measure GRB spectra from 2 keV up to 1 MeV. The POET 
spacecraft provides a zenith-pointed platform for maximizing 
the exposure to deep space and spacecraft rotation provides a 
means of effectively dealing with systematics in the polariza- 
tion response. POET provides sufficient sensitivity and sky 
coverage to detect up to 200 GRBs in a two-year mission. 

LEP and GRAPE determine polarization by measuring the 
number of events versus the event azimuth angle (EAA) as 
projected onto the sky. This is referred to as a modulation 
profile and represents a measure of the polarization magni- 
tude and direction of polarization for the incident beam. De- 
pending on the type of polarimeter, the EAA is either the di- 
rection of the ejected photoelectron (LEP) or the direction of 
the scattered photon (GRAPE). The response of a polarimeter 
to 100% polarized photons can be quantified in terms of the 
modulation factor, /i, which is given by: 

r -c 

^max ^min ^ -, ^ 

Where Cmax and Cmin are the maximum and minimum of the 
modulation profile, respectively. The polarization fraction (11) 
of the incident flux is obtained by dividing the measured mod- 
ulation by that expected for 100% polarized flux. The polar- 
ization angle ($0) corresponds either to the maximum of the 
modulation profile (LEP) or the minimum of the modulation 
profile (GRAPE). To extract these parameters from the data, 
the modulation histograms are fit to the functional form: 

C($) = A+Bcos2($-$o) (2) 
The sensitivity of a polarimeter is defined in terms of the 
minimum detectable polarization (MDP), which refers to the 
minimum level of polarization that is detectable with a given 
observation (or, equivalently, the apparent polarization arising 
from statistical fluctuations in unpolarized data). The precise 
value of the MDP will depend on the source parameters (flu- 
ence, spectrum, etc.) and the polarimeter characteristics. At 
the 99% confidence level, the MDP can be expressed as. 



MDP = 



4.29 Rs + Rh 



fiR, 



(3) 
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TABLE 1 
Instrument Parameters 



GRAPE LEP 



Polarimetry 60-500 keV 2-15 keV 

Detectors BGO/plastic scintillator (62) Ne : CO2 : CH3NO2 Gas (8) 
Spectroscopy 15 keV - 1 IVleV 2-15 keV 

Detectors Nal(TI) scintillator (2) as above 

Field-of-View ±60" ±44° 



where Rs is the observed source strength (cts s~ ), Rb is the 
total observed background rate (cts s"'), and t is the observ- 
ing time (s). The ultimate sensitivity, however, may not be 
limited by statistics but by systematic errors created by false 
modulations that arise from azimuthal asymmetries in the in- 
strument. 

At energies from ^ 50 keV up to several MeV, photon in- 
teractions are dominated by Compton scattering. The opera- 
tional concept for GRAPE is based on the fact that, in Comp- 
ton scattering, photons are preferentially scattered at a right 
angl e to the incide nt electric fi eld vector (the pola rization vec- 
tor) (iBloser et al.l l2008. 2006; Jason et al.ll2005l) . If the inci- 
dent beam of photons is polarized, the azimuthal distribution 
of scattered photons will be asymmetric. The direction of the 
polarization vector is defined by the minimum of the scatter 
angle distribution. The GRAPE performance characteristics 
are shown in Table [T] The design of the GRAPE instrument 
is very modular, with 62 independent polarimeter modules 
and 2 spectroscopy modules. Each polarimeter module incor- 
porates an array of optically independent 5x5x50 mm^ non- 
hygroscopic scintillator elements aligned with and optically 
coupled to the 8x8 scintillation light sensors of a 64-channel 
MAPMT. Two types of scintillators are employed. Low-Z 
plastic scintillator is used as an effective medium for Compton 
scattering. High-Z inorganic scintillator (Bismuth Germanate, 
BGO) is used as a calorimeter, for absorbing the full energy of 
the scattered photon. The arrangement of scintillator elements 
within a module has 28 BGO calorimeter elements surround- 
ing 32 plastic scintillator scattering elements. Valid polarime- 
ter events are those in which a photon Compton scatters in one 
of the plastic elements and is subsequently absorbed in one of 
the BGO elements. These events can be identified as a coin- 
cident detection between one plastic scintillator element and 
one BGO calorimeter element. The azimuthal scatter angle is 
determined for each valid event by the relative locations of hit 
scintillator elements. It is not necessary to know where within 
each element the interaction takes place (e.g., the depth of in- 
teraction). It is sufficient to know only the lateral location of 
each element to generate a histogram of photon scatter angles. 

At energies below ~ 50 keV, the most sensitive technique 
for broadband polarimetry is the photoelectric effect. The 
LEP measures the polarization of incident photons with the 
innovative operati on of a Time Projection Chamber (TPC) 
dBlack et"aLll2007h . The LEP polarimeter enclosure consists 
of four dual-readout detector modules each with an isolated 
gas volume contained by a Be X-ray window. Each detector 
module contains two 6 x 12 x 24 cnP (LxWxH) TPCs that 
share a single X-ray transparent drift electrode. Each TPC is 
comprised of a micropattern proportional counter, consisting 
of a shared drift electrode and a high-field gas electron mul- 
tiplier (GEM) positioned 1 mm from a strip readout plane. 
When an X-ray is absorbed in the gas between the drift elec- 




Ep°°=[keV] 



Fig. 1. — Contour of MDPs (i.e., the minimum values of polarization 
that is detectable at the 99% confidence level) of GRAPE (thick lines) and 
LEP (thin lines) for variable i?p, obs and F and for fixed a = —0.2, /3 = 1 .2, 7" = 
20 sec and the incident angle 30 degree, where -Ep. obs is the observed photon 
energy at the spectral peak, F is the time-averaged flux in 2-400 keV, a and 
l3 are the lower and higher indices of the F,^ spectrum, respectively, and T is 
the duration of the prompt emission. The combination of GRAPE and LEP 
enable us to measure polarizations with reasonable sensitivity in very wide 
energy range. 

trode and the GEM, a photoelectron is ejected in a preferential 
direction with a cos^ $ distribution, where $ is the azimuthal 
angle measured from the X-ray polarization vector As the 
photoelectron travels through the gas it creates a path of ion- 
ization that drifts in a moderate, uniform field to the GEM 
where an avalanche occurs. The charge finally drifts to the 
strip detector where it is read out. 

To estimate realistic MDP values for GRBs detected by 
GRAPE and LEP, we perform an analytical calculation for 
LEP and a Monte Carlo simulation for GRAPE using the cur- 
rent instrument configuration (Table[T]i. The input spectrum in 
the calculation and the simulation is a typical GRB spectrum 
which can be described as a smoothly broken power-law spec- 
tra characterized by photon energy at the i^/v spectral peak, 
Ep^obs, and lowe r and higher indice s of the /v spectrum, a and 
/?, respectively (iBand et al.lll993h . (We tr eat the spec t ral in- 
dices of the specific energy flux /v, while landetaD dHH 
define and /Jb as the indices of the photon number flux, 
i.e., a = -(q!b+ 1) and (3 = -((3b + 1), since we will calculate 
the net polarizations by using specific energy fluxes (equation 
(|7]l).) The various /Syj obs and time-averaged flux in 2-400 keV, 
F, are investigated with fixed a = -0.2, /3 = 1.2, and a burst 
duration of T = 20 s. We also assume the incident angle of 
bursts to be 30 degree off-axis. We interpret simulated events 
with n > MDP as 'Il-measurable events'. Figure 1 shows the 
contour of the MDP values in the /Sp obs plane for GRAPE 
and LEP. As can be seen in the figure, with the combination 
of LEP and GRAPE, it is possible to measure the polarization 
of GRBs with Epohs ranging from a few keV to MeV with 
reasonable sensitivity. 

3. THEORETICAL MODELS 

We calculate the linear polarization for instantaneous emis- 
sion from a thin spherical shell moving radially outward with 
a bulk Lorentz factor 7^1 and an opening angle 9j. The 
comoving-frame emissivity has the functional form of j'^, = 
AQf{v')5{t' -tQ)5{r' - t'q), where Ao is the normalization which 
may depend on direction in the comoving frame and other 
physical quantities of the shell and f{v') represents the spec- 
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tral shape. A prime represents the physical quantities in the 
comoving frame. The delta functions describe the instanta- 
neous emission at f = fo and r = tq. The normalization, Ao, 
has units of erg cm"^ str"' Hz~' . Using the spherical coordi- 
nate system {r,9,(p) in the lab frame, wh ere = is the line 
of sight, we obtain the spectral fluence (jGranot et al.1ll999l: 
IWoods & LoeblfT999H!oka"& Nakamurall2001l) : 



/j/ = ^^ j dcj) J d(cos0)rl- 



(4) 



72( 1-/3 cos 6')2' 

where z and di are the redshift and the luminosity distance of 
the source, respectively, and v' = {\ +z)v^{ \ - (3 cos 6). The 
integration is performed within the jet cone, so that it de- 
pends on the viewing angle 6*,,, i.e., the angle between the 
jet axis and the line of sight. The corresponding Stokes 
parameters of the local emission (i.e., the emission from a 
given point on the shell) are given by j'^ = j'l, Hq cos(2x') and 
j'v' - f''^() siii(2x')> where H'q and %' are the polarization de- 
gree and position angle of the local emission measured in the 
comoving frame, respectively. The Stokes parameters of the 
emission from the whole shell can be obtained by integrating 
those of the local emission similarly to the intensity 

2- \_ 1+^ /"^^ /■^..„.m„2 Aof(iy') ^^/cos(2x) 



(1-/3 cos 61)2"" \ sin(2x) 
(5) 

The polarization degree is Lorentz invariant, i.e., IIq = Hq. 
The position angle x is calculated by taking account of the 
Lorentz transformation of the electromagnetic waves, and it 
is measured from a fixed direction, which we choose to be the 
direction from the line of sight to the jet axis. Then by cal- 
culating {I,Q,U} = dv{Iy,Q^,Uy}, we obtain the time- 
averaged linear polarization in the given wavebands {v\ , 1/2]' 



n = 



(6) 



We consider synchrotron and Compton drag (CD) mecha- 
nisms for the GRB prompt emission. In the synchrotron case, 
the magnetic field consists of a globally ordered field, Bord, 
and small-scale random field, Bmd, i.e., B = Bord + Bmd- The 
field Bord may originate from the central engine, while Bmd 
may be produced in the emission region itself. Here we con- 
sider two extreme cases; synchrotron model with an ordered 
field (SO), in which B^.^ ;» (^rnd)' ^ synchrotron model 
with a random field (SR), in which ^ ^ (Bmd)- For the SO 
model, in particular, we assume a toroidal magnetic field. In 
the following sub-sections, we describe Ao, /{i''), Hq, and x 
as functions of (6',(/)) for each model, and calculate the linear 
polarization for given parameters 7, dj, 9y, and z- 

3.1. SO model: synchrotron with ordered field 

The prompt emission of GRBs could be explained by syn- 
chrotron emission from accelerated electrons that have a non- 
thermal energy spectra by some dissipation process within the 
jet, e.g, internal shocks. Synchrotron emission from the rel- 
ativistically moving shell within a globally ordered magnetic 
field results in a net obs erved linear polarization, reflecting 
the direction of th e field (iLyutikov et al.ll2003t iGranotI 120031 
lNakaretal.ll2003h . Let us assume that the jet is permeated by 
a toroidal field. This is a Ukely configuration if a magnetic 
field is advected by the jet with a constant speed from th e 
central engine (e.g.. Spruit et al.ll200Tl:lFendt & Ouved' 2004 . 

A general formula for calculating the observed linear po- 
larization for synchrotron emission from a uniform jet, in 



which the electrons have a single power-law energy spec- 
trum and an isotropic pitch angle distrib ution and the mag- 
netic field is ordered globally, is derived bv lGranotI (12003 ) and 
iGranot &~Ta vlor (2005). Here we adopt their formulation and 
extend it for the electrons having a broken power-law energy 
spectrum in order to repro duce the typical observed spectra 
of GRBs (iBand et al.lll993h . We adopt the following form for 
the radiation spectrum: /{ly') = fix) where x=v' /v[^ and 



, X ioxx< j3-a 

JW - <; jf-/3(^_Q,)/3-ae"-/3 forjc > f3-a. 



(7) 



where ly^, a, and (3 are the break frequency and low-energy 
and high-energy spectral indices of the comoving spectrum, 
respectively. " If we assume that the energy spectrum of 
the electrons and the strength of the magnetic field are uni- 
form in the emitting shell, then we may write Aq = (sin6'g)"''"', 
where 9g is the angle between the direction of the emit- 
ted radiation and the loca l direction of the magnetic field 
(iRybicki & Li ghtmanll 19791) . The local polarization degree is 
given by: 



Ho 



^ J (a+l)/(a+f) forji:<^- 
" ~ \ (/3+l)/(/3+|) forx>/3- 



(8) 



For a globally ordered magnetic field, the Faraday depo- 
lariz ation effect may be stro ng within the emitting region 
(e^ lToma et al . 2008; Matsu miva & lo ka 2003; Sagi v et alj 
2004), but we neglect it here for simplicity. By using a new 
variable y = (j0) , we obtain (see Appendix lA.lt : 



l-yX 4y (fl- cos 0)2 
l+y J (1+^)2 l+fl2_2flcos(; 



1/2 



X = (/)+arctan 



l-y sm4> 
1 +y a — cosc^ 



(9) 



(10) 



where a = 6/9y. Then the formulation of the net polarization 
degree in the observed frequency region {v\ , vi} becomes: 



Al+qr\j dy 



f{\+qfyi dy fA0Cv) 



where x = {l +z)i^(l +y)/2ji/l^, and 



(11) 



(12) 



A0(y) = 



cos 



(r-i)y;+y 



for ^ > 1 andy < (l-qfyj, 
forq < 1 andy < (1-^)23;^-, 

otherwise. 



(13) 

The polarization degree, 11, in the waveband [i^i , 1^2] can be 
calculated if the geometrical parameters, yj,q, the spectral pa- 
rameters, 7^0, ct, (3, and the redshift, z, are given. 

' ' In our model the radiation spectrum Q is thought to be produced by the 
broken power-law energy spectrum of electrons: N{-fe) oc 7,,''' for -je < 70 
and N(-fe) oc 7^ ''- for 7^ > 70, where a = {pi — l)/2 and /3 = (p2 - l)/2. 
This formulation also includes the case of pi < 1/3, in which a = —1/3, 
Ao oc (sine'g f/^, and Hq^" = 1/2 for a: < <Granolll2003h . 
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Fig. 2. — Linear polarization degrees in the 60—500 keV band as a function 
of q = 6y/6j, where 9y is the viewing angle of the observer and Bj is the 
jet opening angle, for several values of Vj = (-fOj)^, calculated in the SO 
model (synchi'otron model with globally ordered magnetic field). The other 
parameters are -jUq = 350 keV, a = —0.2, /3 = 1 .2, and z = 1 . 

Figure |2] shows the polarization degree in the 60-500 keV 
band as a function of q for several values of yj. The other pa- 
rameters are ji^l^ = 350 ke V, a = -0.2, f3 = 1 .2, and z = 1 . The 
polarization degree is negligible for q^Q, because in this case 
the local polarization vectors are axisymmetric around the line 
of sight, i.e., x = (see Appendix lA.lb . and the local polariza- 
tions are canceled out. For yj > 1, a high level of polarization 

is obtained fory^.''^^ < ^ < 1 (i.e., 7"' < 6',, < Oj). In this case, 

only a fraction of the emitting shell (i.e., 9 < 7"') is bright 
because of the relativistic beaming effect, and the direction of 
the magnetic field is quite ordered in the bright region. The 
contribution of the emission from high latitude, 9 > 7"', is 
negligible, especially foryy > 100, so that the net polarization 
degree is determined only by the emission from the bright re- 
gion with 6 < 7~' and then it is nearly constant. The results 
of our calculations for the case of a = and yj > 100 are con- 
sistent with the results of Granot (2003|) and iLvutikov et al] 
(l2003h . Foryj < 1, a high level of polarization is obtained for 

^ ~ 1 +y/^^ (i.e., 9y 9j + j~^). In this case, the bright region 
on the emitting shell is small, also. 

The polarization is higher for softer spectra (i.e., larger a 
and f3). For example, for yj = 100, -fi^Q = 350 keV, and z = 1, 
the polarization degree at the plateau for ^ < 1 is ~ 0.28 for 
a = -0.5 and (3 = 0.9, while it is ~ 0.52 for a = 0.4 and /3 = 1 .8. 
This is caused mainly by the dependence of the synchrotron 
polarization on the spectral indices (equation |8]l. The maxi- 
mum polarization degree obtained in the SO model is 2± 0.8 
foryj > 0.01, a < 0.4, and /3 < 1.8. 

3.2. SR model: synchrotron with random field 

If the magnetic field is produced at the shock itself within 
the jet, the directions of the field would be random on a scale 
as small as the plasma ski n depth dOruzin ov & Waxman 1999; 
iMedvedev & Loebll 19991) . It is quite plausible that the direc- 
tions of the field are not completely random, but have sym- 
metry around the direction normal to the shock. The less 
isotropic the magnetic field directions behind the shock, the 
higher the local polarization. We consider the extreme case 
in which the field is random strictly within the plane of the 
shock. In this model, the directions of the local polarization 
vectors on the shell are axisymmetric around the fine of sight 



Fig. 3. — Same as Figure[2] but in the SR model (synchrotron model with 
small-scale random magnetic field). 

(see below), so that no net polarization remains if the visible 
region, 6 < 7"', is wholly within the jet cone. However, if 
the observer views the jet from an off-axis angle and the sym- 
metry is broken a high level of polarization r emains (IWaxmanI 
l2003tlSarilfT99llGhisellini & Lazzatil[T999l) . 

Similarly to the SO model, we adopt the broken power- 
law form of the spectrum; /(z^') = fix), where x = v' /v'q 
and f{x) is given by equation d?]). We assume that the en- 
ergy distribution of the electrons and the strength of the mag- 
netic field are uniform in the emitting shell. The local Stokes 
parameters are given by averaging them with respect to the 
magnetic field directions within the shock plane (see Ap- 
pendix lA!2] i. Thus we may write Aq = ((sin^?^)""^'}, where (} 
represents the average. The local polarization degree is given 
by Ho = no^"((sin6l^)"+i cos(20^))/((sin6l^)"+i), where: 



((sin0^r+i) = i r di^' 
I" Jo 



Ay 



H+y)- 



1 I 
■ cos 1J 



1 (a+I)/2 



(14) 



((sin0i,)"+icos(20^)) = 1/^^77' 



4v 2 / 

COS rj 



(I+v> 



(a-I)/2 



sin^ 7/ 



2 / 

cos 77 



(15) 

The local polarization position angle measured in the lab 
frame is given by x = 0^ therefore we obtain the formula- 
tion for the net polarization in the observed frequency region 



n = 



X ((sin6i^)"+i cos(20Jj)) sin(2A0(y)) 



^f(x){(sm9'sr^') 2A0(y) 



1-1 



(16) 



syn 



where q = e,./9j, yj = {-i9jf, x = (1 +z)i/(l +y)/27i/^, and 11^ 
and A0(3') are given by equations (O and ( fT3b . respectively. 

Figure [3] shows the polarization degree in the 60-500 keV 
band as a function of q for several values of yj. The other 
parameters are 71/,, = 350 keV, a = -0.2,(3 = 1.2, and z = 1. 
The results of our ca lculations for t he ca se of a = are con - 
sistent with those of lGranotI ( l2003h and iNakar etall (l2003h . 

-1/2 

A high level of polarization is obtained for ^ ~ 1 +3;^ ' (i.e., 

6*,, ^ 9j + 'y~^) for each value of yj. Since the local polarization 
vectors are axisymmetric around the line of sight, the local 
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polarizations are canceled out if the line of sight is within the 
jet cone. If the jet is observed from an off-axis angle, the net 
polarization remains. The local polarization degree is high- 
est for emission where = 7"' , so that the net polarization 
has a maximum value. The maximum 11 is higher for smaller 
jj, because the contribution of the emission from high lati- 
tude points (6 > 7~'), with a low level of local polarization, is 
smaller 

Similarly to the SO model, the polarization is higher for 
softer spectra, mainly because of the dependence of the local 
polarization degree on frequency (equation [8]l. For example, 
for yj=l, ^v'q = 350 keV, and z = 1 , the maximum polarization 
is ~ 0.32 for a = -0.5 and (3 = 0.9, while it is ~ 0.49 for 
a = 0.4 and /3= 1.8. For > 0.01, a < 0.4, and /3 < 1.8, the 
maximum polarization degree in the SR model is ~ 0.8. 

3.3. CD model: Compton drag model 

The prompt emission from GRBs could be produced by 
bulk inverse Comptonizat i on of soft photo ns from the rel- 
ativistic jet (iLazzati et al I l2004t fe ichler & LevinsonI 120031: 
iLevinson & Eichleil2004tlShaviv & Dara995i r The local po- 
larization position angles are symmetric around the line of 
sight, similarly to the SR model. Therefore this model also 
requires an off-axis observation of the jet to achieve a high 
level of polarization. However, the CD model is different 
from the SR model in the fact that the CD model can in princi- 
ple achieve 11 ^ 1 under the most optimistic geometric config- 
urations, whereas the maximum nis~(/3+l)/(/3+|)^0.8 
in the SR model. 

We assume that the seed radiation is unpolarized and has 
a non-thermal, isotropic spectrum, and the scattered radiation 
has the broken power-law spectrum f{v') = fix), where x = 
v' I v'f^ and /(x) is given by equation (|7]). If the intensity of the 
seed radiation and the electron number density of the shell are 
assumed to be uniform then we may write Aq = (l + cos^6'')/2, 
and Hn = (1 -cos^ 6l')/(l + cos^ Q') (iRvbicki & Lightmanll979t: 
iBegelman & Sikorall987h . The polarization vectors in the co- 
moving frame are perpendicular to both incident and scatter- 
ing directions of photons, so that we obtain %=</)+ ^ in the 
lab frame. Therefore we achieve the formulation for the net 
linear polarization in the observed frequency region \v\ , V2\: 



n = 



/- dvl'^'^'-^' j^fix)^ sin(2A0W) 



(I+9)'yj dy 



(1+V)2 



'(1+y): 



(17) 



where q = 6v/9j, yj = (j9jf, x = (1 +z)j/(l +y)/2-/iyl^, and 
Aip(y) is given by equation (fT3] ). 

Figure |4] shows the polarization degree in the 60-500 keV 
band as a function of q for several values of yj. The other pa- 
rameters are 71^0 = 350 keV, a = -0.2, /? = 1 .2, and z = 1 . The 
results of our calculations for the case of a = /3 are consistent 
with those of ILazzati et al.l (|2004|) . The results are similar to 
those of the SR model, but the polarization degree is higher 
than in the SR model. 

The polarization is higher for softer spectra, although the 
local polarization degree is not dependent on the frequency in 
this model. For instance, for yj = 1, 71^0 = 350 keV, and z = 1, 
the maximum polarization is ~ 0.66 for a = -0.5 and f3 = 0.9, 
while it is ~ 0.7 1 for a = 0.4 and f3= 1.8, but the variation is 
smaller than for the synchrotron models (see § l3.1l and l3.2l l. 
This variation is caused by the kinematic effect. The local 
polarization degree is a maximum for 6 = 7"' (i.e., 6' = 7r/2). 




Fig. 4. — Same as Figure[2] but in the CD model (Compton drag model). 

Thus the net polarization is higher when the contribution of 
the emission from higher latitude with 9 > 7"' is smaller. The 
high latitude emission is dimmer as the radiation spectrum 
is softer. Therefore the net polarization is higher when the 
spectrum is softer. This effect also arises in the SO and SR 
models, although in those models the intrinsic dependence of 
polarization on the spectrum (equation[8]) is rather strong (see 
§[TI]and[321). For yj > 0.01, a < 0.4, and /5 < 1.8, the max- 
imum polarization degree for the CD model is ~ 1 .0. 

4. STATISTICAL PROPERTIES 

In this section we show the results of our Monte Carlo sim- 
ulation of the GRB prompt emission polarization. First, in 
§ 14.11 we give the values of the model parameters so that the 
observed fluences and peak energies of simulated bursts are 
consistent with the data obtained with the HETE-2 satellite. 
In § 14.21 we examine the properties of the polarization distri- 
bution of bursts detectable by the POET satellite, regardless 
of instrument MDR Next, in § 14.31 we show the distribution 
of polarizations that can be measured by POET, and discuss 
how we may constrain the emission models. 

4.1. Model parameters 

We performed Monte Carlo simulations to obtain the dis- 
tribution of the observed spectral energies and fluences in 
the three emission models. Such simulations have been de- 
veloped to discuss the empirical correlation between spectral 
peak energies in the cosmological rest frame and isotropic 7- 
ray energies among GRBs and X-ray flashes in several mod- 
els of geometrical structure of GRB je ts (Zhana et al. 2004t 
'Yamazaki et al."2004'; 'Dai & Zhang 2005; Toma et al. 20051 
Donaghy 2006). We generated 10,000 GRB jets with Lorentz 
factor, 7, and opening angle, Oj, and a random viewing an- 
gle for each jet according to the probability distribution of 
sin6yd9i,d(j) with 6*,, < 0.22 rad. For each burst generated 
we calculate the i^I^ spectrum to obtain the spectral peak en- 
ergy, /ip.obs, and the fluence, /, in the 2-400 keV range by 
using equation (|4]l. Since obs's and /'s calculated for each 
q = 9y/6j in the three models are different only by factors less 

We confirmed that the bursts with 9v > 0.22 rad in our simulation 
are not detected by HETE-2 or POET with the parameters we adopt in this 
paper We can therefore discuss the distribution of several quantities of the 
detectable bursts and the event rate ratio of bursts for which polarizations 
can be measured to the detectable bursts without considering the bursts with 
0, > 0.22 rad. 
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Fig. 5. — The Epabs-F diagram calculated in our Monte Carlo simu- 
lation. The simulated events that can be detected by WXM on HETE-2 are 
represented by dot s. They are co mpared with the HETE-2 data (points with 
errorbars) (Sakamoto et aU2005D . 

than 2, Ep obs's and /'s of the simulated bursts may be calcu- 
lated using just one model, for which we chose the CD model. 

The distributions of 7 and 6j for GRB jets are highly uncer- 
tain. We make a simple assumption for the distribution and in 
§ |4.3| we perform some simulations for different assumptions. 
We fix 7 = 100. We assume the distribution of 9j as 



f(6j)d6j « 



0fd0j 



for 0.001 <6I^ <0.02 
for0.02<6l, <0.2, 



(18) 



where qi = 0.5 and ^2 = -2.0. The value of ^2 = -2 is in- 
ferred from the observations of the steep ening breaks (i.e. , 
jet breaks) of some optical afterglows dFrail et aTl 120011 : 
IZeh et al.ll2006h and from analysis of BAT SE data using some 
empirical relations (lYonetoku et al.ll2005h . The re are several 
suggestion s of even ts with very small 9j (e.g., ISchadv et alj 
12007 : R acusin et al.ir2008) . although the value of qi is highly 
uncertain. The spectral parameters rlAo,-fi'Q,a, and /? are as- 
sumed as follows. The first two parameters are given so that 
the rest-frame spectral peak energies and isotropic 7-ray en- 
ergies calculated for a jet viewed with 9y = are consistent 
with those of typical GRBs. Such an on-axis emission has ap- 
proximately Ep = and E-^^o = 167r^rgAo7i'o- The param- 
eters TqAq and ^v'q are given through the empirical relations 
E,,J ln = 10^'6 erg and E„ = 8()6.(£isn /10^' erg)i/2 keV 
(e.g.. iFrail e"tani2001l: lAmati et alJl2002h . We assume that 
the coefficients £1 and £ 2 obey the log-normal distribution 
jloka & Nakamurall2002h with averages of 1 and logarithmic 
variances of 0.3 and 0.15, respectively. The last two parame- 
ters are fixed by a = -0.2 a nd 0=1.2, which are typical values 
for G RB prompt emission (iPreece et alJ2000l:ISakamoto et alJ 
l2005h . The distribution of the source redshift, z, is assumed to 
be in proportional to the cosmic star f ormation rate. We adopt 
the model SF2 in Porciani & Madaul (|2001), i.e., the comov- 
ing GRB rate density is assumed to be proportional to 



exp(3.4z) V»m(1+z)^ + 0a 
exp(3.4z) + 22 (1 + 2)^/2 



standard cosmological parameters 
"\rjM = 0.3,and Ha =0-7 ■ 



(19) 



of //n = 



We take the 

70 km s"' Mpc~ , s = 

Figure|5]shows the results of fip obs and time-averaged flux, 
F. The time-averaged flux is calculated hy F = I/T, where 
T is the duration of a burst. We fix T = 20 s, which is a typ- 
ical value for long GRBs (e.g.. ISakamoto et allllOOl) . We 



Fig. 6.— E„ 

obs ~ n diagrams for the simulated events that can be detected 
by GRAPE in the SO (red open circles), SR (green filled circles), and CD 
(blue plus signs) models. The adopted parameters are as follows. The fixed 
parameters are 7 = 100, q[ = 0.5,^2 = -2.0, a = -0.2, /3 = 1.2, and T = 20 s. 
The distribution of the source redshift z is assumed to be in proportional to the 
cosmic star formation rate. The parameters r^Ag and ■yUg are distributed so 
that the simulated -Ep. obs diagram is consistent with the observed data (see 
Figure|5). See text for the cases of the spectral indices distributed realistically, 
for -0.5 < Q <0.4 andO.9 </3 < 1.8. 

show only the simulated bursts that have fluxes above the de- 
tectable limit of the HETE-2 satellite. They are consistent 
with the data obtained by HETE-2 (Sakamoto et al. 2005). 
The scatter of the simulated bursts is due both to the scat- 
ter of the assumed jet parameters and t h 
dYamazaki et al.ll20o4lDonaghvll2006l) . 



arameters and t he viewing angle effect 

' " ■ 13 



4.2. Properties of polarization distribution 

We calculated the linear polarization, 11, by using equations 
( fTTT i. ( fTSI l. and ([TtT i to obtain the polarization distribution of 
the simulated bursts that can be detected by GRAPE and LEP. 
The detection limits of GRAPE and LEP are set to be the 
MDP contours of 1.0. (see Figure [T]! Figures |6] and [T] 
show the Ep ohs - H diagrams of all the simulated bursts that 
can be detected by GRAPE and LEP, respectively, in the SO 
(red open circles), SR (green filled circles), and CD (blue plus 
signs) models. In the SO model, most of the detectable bursts 
have 0.3 < n < 0.5 in the GRAPE band (60-500 keV), while 
they have 0.2 < H < 0.3 in the LEP band (2-15 keV). In the 
SR and CD models, most of the detectable bursts have 11 < 0. 1 
in both GRAPE and LEP bands. The events with H > 0.1 are 
distributed uniformly with 11 < 0.4 and 11 < 0.9 for the SR 
and CD models, respectively. 

This result can be roughly explained by the polarizations 
calculated as functions of y ; and q = (?,,/ dj for = 350 keV 
and z = 1 (see Figures |2] [3] Hll and the distribution of dj and q 
for the detectable bursts in this simulation, shown in Figure [8] 
The detectable events are dominated by the events with ^ < 1 , 
since events with < 1 are much brighter than those with 
q > I because of the relativistic beaming effect. For events 
with q > 1, narrower jets are easier to detect since they have 
intrinsically higher emissivities by our assumption. Most of 
the detectable events have ^ < 1 and dj > 0.02 (i.e., yj > 4). 

" lYamazaki et alj 120041) showed a deviation from the Amati relation 

1/2 

(Ep oc E-J^ ) for Ep < 10 keV in the uniform jet model, but the Ep.obs — E 
diagram we derive is still consistent with the observed dataset. 

The detection limits of GRAPE and LEP for signal-to-noise ratio > 5 
are similar but not identical to the MDP contours of 1.0. Thus our setting for 
the detection limits is just for simplicity. 
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Fig. 7. — Same as Figure|7] but for LEP. 
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Fig. 8. — Distribution of 9j and q = dyjdj of the detectable bursts by 
GRAPE in the model described by Figure|6] 

For these events the SO model gives 0.3 < 11 < 0.5 in most 
cases, while the SR and CD models give 11 < 0.1, for the 
GRAPE band as shown in Figures |2][3] and|4] The remaining 
detectable events mainly have ^ > 1 and Bj > 0.005 (i.e., yj > 
0.25). These events have 11 < 0.6 in the SO model, 11 < 0.5 in 
the SR model, and H < 0.9 in the CD model, for the GRAPE 
band as shown in Figures |2][3] and|4] The results for the LEP 
band can be explained similarly. 

In all the three models, the results show 11(60- 500 keV) > 
n(2- 15 keV) for almost all the detectable bursts with 11 > 
0.1. This is due to the fact that typically the contribution of 
the high-energy photons with spectral index (3 is larger in the 
GRAPE band than in the LEP band. The emission with softer 
spectrum has higher polarization because of the intrinsic prop- 
erty of the synchrotron polarization (equation [H]) for the SO 
and SR models and the kinematic effect for the CD model 
(see § I3.3l l, respectively. 

In the SO model, the polarization of GRBs with q < I is 
higher for lower /Syj.obs for the GRAPE band. This is because 
the contribution from high-energy photons, with energy spec- 
tral index /3, is larger In the SR and CD models, the higher 
n GRBs can be obtained for smaller 9j. The maximum 11 is 
obtained for Oj ~ 0.002. 

4.3. Cumulative distribution of measurable polarizations 



We obtain the distribution of polarization that can be mea- 
sured, by using the MDP values we derived for a = -0.2, 
/3 = 1.2, and T = 20 s (see § We interpret the simulated 
events with 11 > MDP as 'Il-measurable events'. Figure |9] 
shows the cumulative distribution of 11 that can be measured 
by GRAPE and LEP in the SO, SR, and CD models. We 
have set the number of detectable events Nj = 200. Since 
the polarization in the LEP band is lower than in the GRAPE 
band for almost all the cases as discussed in § 14.21 the num- 
ber of events for which polarization can be measured by LEP 
is smaller than for GRAPE. In the SO model, the number of 
Il-measurable bursts is N„, > 60, and the cumulative distribu- 
tion of measurable 11 varies rapidly at 0.3 < 11 < 0.4 for the 
GRAPE band. In the SR model. A',,, < 10, and the maximum 
polarization is Ilmax < 0.4. In the CD model. A',,, < 30, and 
n„,ax < 0.8. 

To investigate general properties of the cumulative distri- 
bution that do not depend on the model parameters, we per- 
formed simulations for other values of 7, qi, and q2, the 
Lorentz factor of the jets and the power-law indices of the 
distribution of the opening angles of the jets, respectively. We 
refer to the parameters adopted for the above simulation as 
'typical' parameters. We now consider a range of parameters: 
7 > 100, ^1 > 0.5, and ^2 > ~3 0. which are quite reasonabl e 
for GRBs (e.g..' Lithwick & Sarill2001l:lYonetoku et al.l |2005l). 
Within these parameter ranges we obtain the lower (upper) 
limit of N„,/Nd for the SO model (the SR/CD models). 

Figure [To] shows the results for 7 = 300 and the same 'typi- 
cal' values for the other parameters. The number A^„, is larger 
in the SO model and smaller in the SR and CD models than 
the case of 7 = 100. As 7 is larger, the beaming effect is 
stronger and the ratio of the bursts with ^ < 1 for detectable 
bursts is larger. Thus the number of bursts with a high degree 
of polarization is larger in the SO model and smaller in the 
SR and CD models. Figure [TTI shows the results for ^1 = 1.0 
and the same 'typical' values for the other parameters. Since 
the ratio of the number of the bursts with smaller v/ to that of 
detectable bursts is smaller. A',,, is slightly smaller than that for 
the 'typical' parameters in the SR and CD models. Figure [T2l 
shows the results for q2 = -3.0 and the 'typical' values for the 
other parameters. In this case A',,, is slightly larger than that 
for the 'typical' parameters in the SR and CD models. The 
number A',,, in the SO model is similar for Figure |9] [H] and 
Figure[T2]in the GRAPE band. To summarize, for the parame- 
ters 7 > 100, ^1 > 0.5, q2 > -3.0, a = -0.2 and /? = 1 .2, we can 
say thatN,„/Nd > 30% for GRAPE and the cumulative distri- 
bution of measurable 11 varies rapidly from 0.3 < 11 < 0.4 in 
the SO model. For the SR model, N„,/Nd < 5% for GRAPE, 
with a maximum polarization Ilmax < 0.4. For the CD model, 
N,n/Nd < 15% for GRAPE, and n^^ax < 0.8. 

Since the dependence of the polarization degree on the 
spectral indices is relatively large in the SO and SR mod- 
els, we should take account of the distribution of a and f3. 
The observed spectral parameters a and (3 ar e distributed 
roughl y as -0.5 < a < 0.4 and 0.9 < /3 < 1.8 (iPreece et alJ 
2000; Sakamoto et al.ll2005l) . Within these ranges of a and 
(3, the polarization degree for yj > 10, q < I, and 50 < 
Epochs < 10^ keV is 0.2 < n < 0.7 in the SO model. Thus 
the measurable polarizations are clustered at 0.2 < 11 < 0.7. 
The maximum polarization obtained in the SO mode l for 
yj > 0.01, a < 0.4, and /? < 1.8 is ~ 0.8 (see § [TB. In 
this case N,„/Nd will be larger than 30%. In the CD model, 
the result will not be significantly different from the case of 
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fixed a and /?. In the SR model, the polarization degree 
does not exceed those calculated in the CD model, and thus 
Nm/Nd < 15%. The maximum polarization obtained in the 
SR model for yj > 0.01, a < 0.4, and /? < 1.8 is ~ 0.8 (see 
§E3. 

In conclusion, we can constrain the emission mechanism 
of GRBs by using the cumulative distribution obtained by 
GRAPE. If N,n/Nd > 30%, the SR and CD models may be 
ruled out, and in this case if the measured polarizations are 
clustered at 0.2 < 11 < 0.7, the SO model will be favored. If 
Nm/Nd < 15%, the SO model may be ruled out, but we can- 
not distinguish between the SR and CD models with different 
distributions of yj, a, and f3. If several bursts with 11 > 0.8 
are detected, however, the CD model which includes adequate 
number of small yj bursts will be favored. 

5. SUMMARY AND DISCUSSION 

Recently there has been an increasing interest in the mea- 
surement of X-ray and 7-ray polarization, and the observa- 
tional techniques can now achieve significant sensitivity in 
the relevant energy bands. Several polarimetry mission con- 
cepts, such as POET, are being planned. The POET con- 
cept has two polarimeters, GRAPE (60-500 keV) and LEP (2- 
15 keV) both of which have wide fields of view. If launched, 
missions of this type would provide the first definitive de- 
tection of the polarization of GRB prompt emission. This 
would enable the discussion of the statistical properties of 
the polarization degree and polarization spectra, which will 
give us diagnostic information on the emission mechanism 
of GRBs and the nature of the GRB jets that cannot be ob- 
tained from current spectra and lightcurve observations. We 
have performed Monte Carlo simulations of the linear po- 
larization from GRB jets for three major emission models: 
synchrotron model with globally ordered magnetic field (SO 
model), synchrotron model with small-scale random magnetic 
field (SR model), and Compton drag model (CD model). We 
assumed that the physical quantities for the emission of the 
jets are uniform on the emitting surface and that the jets have 
sharp edges. Our jet angle distribution allows the detections 
of GRBs with very small opening angles (i.e., smal ler than 
1 degree) as suggested by several Swift bursts ( Schad y et alJ 
I2007t iRacusin et alj|2008h . We have shown that the POET 
mission or other polarimeters with similar capabilities, i.e., 
broadband spectral capabilities for the determination of Ep ^bs 
and sensitive broadband polarimetric capabilities to minimize 
MDP, can constrain the emission models of GRBs. Further- 
more, these simulations indicate that an increase in the LEP 
effective area would be beneficial to compensate for the lower 
expected polarization at lower energies. 

As shown in Figures |2l [3] and |4] the SR and CD models 
require off-axis observations of the jets to achieve a high level 
of polarization, while the SO model does not. In this sense the 
SR and CD models are categorized a s geometric models, and 
the S O model as an intrinsic model (IWaxmanll2003l : iLazzatil 
12006*). The distribution of observed polarizations obtained by 
our simulations show that the geometric SR/CD models will 
be ruled out if the number ratio of the Il-measurable bursts 
to detected bursts is larger than 30%, and in this case the SO 
model will be favored if the measurable polarizations are clus- 
tered at 0.2 < n < 0.7. If the number ratio is smaller than 
15%, the SO model may be ruled out, but we cannot distin- 
guish between the SR and CD models with different distribu- 
tions of yj = i'^Oj)^, a, and (3, where 7 and 9j are the bulk 
Lorentz factor and the opening angle of the GRB jet, respec- 



tively, and a and f3 are lower and higher indices of the energy 
spectrum. However, if several bursts with 11 > 0.8 are de- 
tected, the CD model which includes an adequate number of 
small yj bursts will be favored. 

If the cumulative distribution of the measurable polariza- 
tions favors the SO model, the globally ordered magnetic field 
would be advected from the central engine. If we understand 
the strength of the magnetic field in the emitting region from 
the luminosity and the spectrum of the emission, we can con- 
strain the strength of the field at the central engine. If the 
geometric SR/CD models are favored from the observations, 
it will be established, independently of the afterglow obser- 
vations, that GRB outflows are not spherical but highly col- 
limated. If the CD model is favored by the observations, we 
may constrain the distribution of the parameter yj = (•yOj)^ of 
GRB jets. The CD model needs a dense optical/UV photon 
field interacting within the relativistic jets jLazzati et al.l2000l : 
lEichler & Levinsonll2003h . 

We have made some simplifications in our simulations, and 
there are some caveats. We have assumed that the jets are 
uniform on the emitting surfaces and have sharp edges. To 
compare the simulations and the observa tions further , more 
sophisticated mo dehng is required (e.g.. IZhang et al.l 120041 : 
iToma et al.ll2005h . 

We have interpreted bursts as a simple combination of 
pulses, without taking account of the temporal variation of the 
Lorentz factor 7 of the jet. If this is accounted for, each pulse 
may have different yj = {'jOj)^ but the same q = 9^/0 j- We 
should then average the polarizati on with respect to fluence 
of ea ch pulse having different yj (iGranotI 120031 : |Nakar et alJ 
|2003h . However, in the SO model, the cumulative distribution 
of measurable 11 will not be changed significantly as long as 
yj > 10, because 11 is clustered into a small range for ^ < 1 
and yj > 10. To average the polarization in the case of ^ > 1, 
the relation between the luminosity and the Lorentz factor for 
each pulse is required to predict the polarization distribution. 

For the SR model we have assumed that the directions of the 
magnetic field are confined within the shock plane. They may 
be more isotropic in reality, in which case the polarization 
degree in the SR model will be reduced. 

In the synchrotron model with a combination of the glob- 
ally ordered magnetic field and the locally random field, 
B = Bord + Brnd, the linear polarization can be calculated by 

n = (gord + 6md)/(/oid + /rnd) « (Ilord + V^md) / ( 1 + ??), whcrC 

{/, 2}oi-d and {/, 2}rnd are the Stokes parameters from the or- 
dered and random fields, respectively. Ilord and Xlmd are de- 
scribed by equations (fTTI) and ( fT6b , and rj = (Bmd/fiord)"^'- 
This model will reduce the number ratio of Il-measurable 
bursts to detected bursts to less than 30% and the clustering 
of measurable polarizations will be at 11 < 0.7. 
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Fig. 9. — The cumulative distribution of 11 that can be measured by GRAPE (left) and LEP (right) in the SO (solid), SR (dashed), and CD (dot-dashed) 
models in which the number of detectable bursts is 200. The adopted parameters are as follows. The fixed parameters are 7 = IQO, tji = 0.5, ^2 = —2.0, a = -0.2, 
/3 = 1.2, and T = 20 s. The distribution of the source redshift z is assumed to be in proportional to the cosmic star formation rate. The parameters r^Ag and ■yUg 
are distributed so that the simulated obs — ^ diagram is consistent with the observed date (see Figure|5). See text for the cases of the spectral indices distributed 
realistically, for -0.5 < a < 0.4 and 0.9 < /3 < 1.8. 
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APPENDIX 

SOME NOTES ON SYNCHROTRON POLARIZATION 

The SO model: synchrotron with ordered field 

We consider the synchrotron radiation from the shell moving radially outward with a bulk Lorentz factor 7^1, and the 
magnetic field in the shell is globally ordered within the plane parallel to the shock plane. If the matter of the shell expands with 
a constant speed, the strength of magnetic field with radial direction scales as while that with transverse direction scales as 
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Fig. 10. — Same as Figure|9] but the Lorentz factor of the jets 7 = 300. 
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Fig. 1 1 . — Same as Figure|9] but the lower power-law index of the 9j distribution ^1 = 1.0 (see equation[T8). 
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Fig. 12. — Same as Figure|9] but the higher power-law index of the dj distribution qi = —3.0 (see equation fTsl . 



. Thus the field advected with the shell is likely to have the direction parallel to the shock plane. 
We set the line of sight (i.e., the direction from the central engine to the earth) in the lab frame to be z axis, and the direction 
of the magnetic field on a given point of the shell, projected onto the plane perpendicular to z axis, to be x axis. The given 
point can be described by spherical coordinates (0, ip). Then the components of the velocity vector of the given point and the unit 
wave vector can be described by the right-handed coordinate system xyz as (3 = { f3 sin d cos ip, f3 sin 6 sin (p,P cos 6) and k = (0,0, 1), 
respectively. The unit wave vector in the comoving frame is 



1 



7(l-/?.k) 



k+/3 



T 



j+ 1 



/3.k-7 



(Al) 



Since the direction of the magnetic field in the comoving frame is perpendicular to the velocity vector of the fluid, B' = 
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(cos B I \/ cos^ Q + sin^ Q cos^ 0, - sin cos kpj \J cos^ Q + sin^ cos^ </?). Then we may calculate cos = B' • k', and we obtain 

sin6l!;«(- I cos^o9 + sin^o9 (A2) 

\ 1 +7^6'^/ 

in the limit 7 ;2> 1 . 

The direction of the polarization vector of the synchrotron radiation is calculated by e' || B' x k'. Then we obtain the direction 
of the polarization vector in the lab frame by 

e = 7(l+^-k')e'-(/3-e') ^-^/3 + 7k'j . (A3) 

The results are = and 

e,. P-cos9 1 
tanxfi = -^=tan(p- ^ . • (A4) 

e.f iSsin 9 sm (f cos ip 

The angle xb is the polarization position angle measured from the axis x (i.e., the direction of the local magnet ic field) . The a bove 
equation can be rewritten as xb ~ (y5 + arctan[(l -7^6'^)cot(/p/(l + 7^6*^)]. This result is consistent with that of iGranod (l2003h . 

Based on the above results, we consider the case that the magnetic field is axisymmetric around the jet and has a toroidal 
configuration. We set the direction from the line of sight to the jet axis to be x axis. Below we will rewrite the above results by 
using the azimuthal angle (f) measured from x axis. In the coordinate system of xyz, the jet axis and the coordinates of a given 
point on the shell are described as J = (sin 6*,, , 0, cos 6',,), and R = (sin 6 cos cj), sinO sin 0, cos 9), respectively. The magnetic field at 
the given point is given byB' = RxJ/|RxJ|. Let the unit vectors of the directions of R and B' projected onto xy plane be f and 
b, and cos 1^9 = b • r. Then we obtain 

2 sin^ (f) 

1 + - 2fl cos (/) 

where a = 9/9y. Equation (|9]l is given by inserting equation dASI l into equation ( IA2b . If we measure the position ang le from the 
X axis, we obtain equation JTOl l, i.e., X = Xb - 0. These results are consistent with those of iGranot & Tavlo^ (l2005h . 

The SR model: synchrotron with random field 

Here we consider that the directions of the magnetic fields are confined within the plane parallel to the shock plane and that 
they are completely random. This field configuration is possible if the field is generated by the shock. In the comoving frame of 
the shell, we set the direction of k' to be axis 3, and set a right-handed coordinate system 123. Let the polar and azimuthal angles 
of B' be 9g and 0^, respectively. In this coordinate system, the Stokes parameters of synchrotron emissivity are given by 

= -j'l,nQ cos(20;j), jlH = -j'l,Ilo sind^'s). (A6) 

Next we set another right-handed coordinate system xyz of which z axis is along the velocity vector of the fluid and xz plane 
includes k'. Then the angle between k' and z axis is 9'. Here the magnetic field B' is confined within xy plane. Let the azimuthal 
angle of B' be ?/, and we obtain the relations between the components of B' in the systems 123 and xyz. 

sin d'g sin (f)'g = cos 9' COST]' , 

sin 9'g cos<j)g = sin r\' , ( A7) 

cos9g = sin 0' COS 77'. 

Then we obtain 

11/2 



sin 61^= l-^y^^^cos^ry'^ 
cos(2</,i,)=-4^pry'-(|^)'cos2 



(A8) 



To obtain the polarization degree of synchrotron radiation from the random field, we average the Stokes parameters 
with respect to ij'. This leads to (f^,) = 0. Then we can calculate the polarization degree by Hq = (f^) / {j'^r) = 
nQ*^"((sin6'g)"'^' cos(20^)) / ((sin6'g)"'^'), and the polarization vector is along axis 1, i.e., the direction perpendicular to k' and 
within the plane including k' and f3. 



